Jamming transitions in amorphous packings of frictionless spheres occur over a 

continuous range of volume fractions 
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We numerically produce fully amorphous assemblies of frictionless spheres in three dimensions 
and study the jamming transition these packings undergo at large volume fractions. We specify four 
protocols yielding a critical value for the jamming volume fraction which is sharply defined in the 
limit of large system size, but is different for each protocol. Thus, we directly establish the existence 
of a continuous range of volume fraction where nonequilibrium jamming transitions occur. However, 
these jamming transitions share the same critical behaviour. Our results suggest that, even in the 
absence of partial crystalline ordering, a unique location of a random close packing does not exist, 
and that volume fraction alone is not sufficient to describe the properties of jammed states. 

PACS numbers: 05.10.-a, 05.20.Jj, 64.70.Pf 



Everyday experience shows that it is not possible to 
compress a disordered assembly of rigid particles beyond 
a maximal packing fraction. The idea of a critical value 
for the volume fraction of this 'random close packing' 
has a long history in science For the idealized case 
of hard, spherical, frictionless particles the idea recently 
emerged that a critical value of the volume fraction, 
named can be appropriately defined in the thermo- 
dynamic limit [3, Q , and that it corresponds to a critical 
point ('point J') with remarkable scaling properties ob- 
served when J is approached from either side 0, 13| ■ ^ 
unique critical packing fraction is also found under quasi- 
static [ij or continuous 0, 0] shearing. Further studies 
have also established the peculiar properties of config- 
urations prepared exactly at with specific, singular 
behaviour of the pair correlation function d, S. This 
leads to the conjecture that the properties of this criti- 
cal point can influence the physical behaviour of dense 
particle systems under various conditions [lol |. 

The identification of a unique critical point for jam- 
ming is surprising because two distinct phenomena occur 
when compressing hard spheres at equilibrium. First, an 
equilibrated system of monodisperse hard spheres is fully 
crystalline above (f) ~ 0.54. Since crystalline states are 
efficiently packed, it is in principle always possible to in- 
crease the volume fraction of a disordered assembly by 
increasing the local order ll|. In the studies mentioned 



above, crystallization is prevented by using size polydis- 
persity which increases the nucleation barrier and can 
efficiently suppress crystallization at equilibrium. More- 
over, these studies employ athermal numerical protocols 
which do not allow the system to reach equilibrium. 

A second relevant phenomenon occurring at equilib- 
rium is the glass transition. Compressing hard spheres 
at thermal equilibrium, it is found that ergodicity is lost 
near a volume fraction much lower than because the 



relaxation time for structural relaxation becomes larger 
than experimental timescales [1^. Necessarily, then, 



the properties of hard sphere glasses should depend on 
the preparation history, as is well-known for molecular 
glasses [l3|. The resulting interplay between glass and 
jamming transitions is currently receiving attention, both 
from theory and simulations. A special density akin to (/),; 
appears in mean-field glasses undergoing a random first 
order transition Il4| . In this language, corresponds to 



the (sharp) appearance of an exponentially large num- 
ber of metastable states, and such states exist along a 
continuous range of densities between 0j and that of the 
'glass close packing' [13, [Hj, where the configurational 
entropy (or complexity) counting metastable states van- 
ishes. However, none of these mean-field concepts is ex- 
pected to remain sharply defined in three dimensions [31 , 
and it is important to assess the validity of this interpre- 
tation in finite dimensions. 

Previous numerical work has shown that different com- 
pression rates produce glasses with a pressure which ap- 
pears to diverge at different densities, at least for fi- 
nite size systems [3, [13, H, Indeed, the analogy 
of the jamming density with the energy of inher- 
ent structures in systems with soft potentials [l3|, sug- 
gests that the jamming densities should depend both on 
compression rates as well as the initial state from which 
the glasses are generated. Recently, a compression pro- 
tocol was specifically devised to prevent any ordering, 
thus allowing to focus directly on the sole influence of 
glassy behaviour on jamming |22l |. Starting from a well- 
equilibrated fluid configuration at a given volume frac- 
tion, 0init, one uses a very large compression rate to 
reach nearly jammed configurations. During these rapid 
compressions, the glass pressure was observed to diverge 
within a finite ran ge o f volume fractions, even in the ther- 
modynamic limit [221, but the jamming transition was 
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not studied. 

Since evidence is mounting in favour of both a unique 
location of the jamming transition with critical properties 
in some papers, or a continuous range of volume fraction 
in some other works, it appears timely to reconcile these 
two lines of research and answer the following impor- 
tant, open questions. Does a continuous range of volume 
fraction for jamming necessarily result from crystalliza- 
tion or demixing? Can one reconcile the results found 
using thermal and athermal protocols? Can the remark- 
able properties of the jamming transition survive if its 
location is not unique? Here, we provide precise answers 
to these questions and close the gap between two sets 
of ideas. We study numerically how amorphous assem- 
blies of hard spheres prepared using the (thermal) tools 
of Ref. [2^ jam at large volume fraction. We then apply 
the alternative (athermal) tools of Ref. 0] to analyze the 
properties of jammed configuration at and near the tran- 
sition. For the specific three-dimensional binary system 
we use, we can directly establish the existence of a con- 
tinuous line of J-points extending at least in the range 
(pj € [0.648,0.662], all sharing similar critical properties 
and exponents. Thus, although point J is not unique, its 
critical properties are. 

We study the jamming transition in a three- 
dimensional system of spherical particles using periodic 
boundary conditions. We use a 50:50 binary mixture 
of spheres with diameter ratio 1.4. Our preparation 
protocol combines the fast compression of hard spheres 
described above [2^, to the athermal cornpression of 
spheres treated as soft spheres as in Ref. [2|. Specifi- 
cally, the particles are first treated as hard spheres and 
equilibrated using Monte-Carlo simulations over a broad 
range of volume fraction up to = 0.596. This allows 
one to obtain disordered configurations representative of 
the metastable fluid of hard spheres (no crystallisation 
or demixing is observed for simulations as long as 10^*^ 
Monte-Carlo timesteps). We then use Monte-Carlo sim- 
ulations to rapidly compress these equilibrated fluid con- 
figurations. During these compressions the system has 
no time to relax and retains a structure very close to the 
initial fluid states. Thus, contrary to previous work em- 
ploying slow compressions [ii, [il, we are certain that 
the jammed configurations we produce contain no more 
order than the original equilibrated fluid configurations. 

When (reduced) pressure is very large, P/{pkBT) ^ 
10^, we stop the compression and switch to the athermal 
procedure of Ref. [1], where spheres are now treated as 
harmonic spheres interacting through a soft pair poten- 
tial, Vsoitinj) = (1 - fij/o-ij)^. Here, a^j = (ct^ -|- aj)/2 
and CTi is the radius of particle i. The compression then 
proceeds in a succession of small instantaneous particle 
inflation followed by energy minimization using conjugate 
gradient • At large volume fraction, the energy cannot 
be minimized to zero anymore, and particles now overlap, 
signalling that the jamming transition has been crossed. 



<^init 


0.3572 0.5397 0.5672 0.5935 


(^J (8000, 0init) 
(^J (1000, 0init) 


0.6481 0.6499 0.6537 0.6616 

0.6466 0.6491 0.6531 0.6616 


a J (8000, 0init) 
a J (lOOO,0init) 


5.2 ■ 10"" 4.9 ■ 10"'' 3.0 • 10"" 2.5 • 10"" 
9.7 ■ 10"" 1.0 • 10"^ 7.6 ■ 10"" 5.8 • 10"" 


2 at (f>j 
212 at f/ij 


6.0017 6.0022 6.0021 6.0023 
2.47 2.50 2.54 2.65 
2.92 2.91 2.92 2.97 


Rattlers 


4.3 % 4.7 % 5.4 % 6.2 % 



TABLE I: Statistics of jammed configurations for difi'erent 
system sizes A'^ and various (\>\a\x- The limit <j)j{N 00, 0init) 
exists, with finite N fiuctuations decaying in good agreement 
with A"^^^, and retains a dependence on 0init' At the tran- 
sition, configurations are all isostatic, 2 ~ 6 (with similar 
distribution width a^), nearly identical local structure, num- 
ber of rattlers, and no increasing demixing between species 
(212 is nearly constant). 



We emphasize that all the details entering this compres- 
sion protocol in principle quantitatively affect the results 
presented below [ist . However, since the most sensitive 
control parameter is (/linit, we shall only use its value to 
distinguish between the different protocols, all other pa- 
rameters being kept fixed. We focus on four values for 
0init listed in Table HI We monitor finite size effects and 
convergence towards the thermodynamic limit by study- 
ing two system sizes, TV = 1000 and 8000, chosen in the 
regime where scaling with N is well understood Q ■ 

In Fig. [T]-a we show the final part (e.g. when particles 
are treated as soft) of a randomly selected compression 
history for a system with N = 10'^. We follow the energy 
density of the system, e = ^^^^Vsoitifij), for in- 
creasing (f>. We have defined configurations having energy 
density less than 10"^^ as unjammed. For low enough 0, 
we obtain only unjammed particle configurations. Upon 
increasing the density, there appears a volume fraction 
above which the energy increases rapidly above zero: this 
is the jamming transition Upon further compression, 
the energy keeps increasing on average, and so does the 
average number of contacts per particle, z. 

This qualitative description becomes more complicated 
when multiple histories are superimposed, see Fig. [T]-b. 
There are several sources of fiuctuations in this plot. A 
first source of fluctuations stems from the possibility for 
the system to undergo some reorganization which makes 
the energy drop suddenly, and makes the e{4>) curve mul- 
tivalued during a single compression run. Thus, no sim- 
ple scaling law for the energy can be detected during 
compressions. The second source of fluctuations arises 
when different realizations of the same protocol (i.e. same 
value of (pinit but starting from independent fluid config- 
urations) are compared. These fluctuations are present 
because N is finite, and should vanish when ^ 00. A 
third source of fluctuations can be seen in Fig. [T] when 
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FIG. 1: (a) The final part of a compression-decompression 
cycle across the jamming transition for A'' — and (^init ~ 
0.3572, showing the practical determination of t/ij for this 
sample, (b) Evolution of the energy density during compres- 
sions of several runs with A'' = 1000. Different curves with the 
same color correspond to independent realizations at constant 
(/>init, while different (phiit are shown with different colors, (c) 
Same as (b) for the evolution of the number of contacts per 
particle, (d) All data in (b) and (c) collapse when (p is elimi- 
nated, showing that the jamming transition is better defined 
by its properties than by its location. 



'/'init is varied, larger 0init typically yielding larger volume 
fraction for the onset of positive energies. 

If we decompress jammed configurations, reorganiza- 
tions are much less likely, and the e{(j}) curve is smooth 
enough that it allows a quantitative determination of 
0j, sec Fig. [T]-a. Of course, (f>j now depends on one 
additional parameter, namely the volume fraction from 
which decompression starts. We arbitrarily decompress 
when e w 10~^ and obtain (j>j by fitting e to a power law 
decay, e ~ (0 — 0./)^ during decompression. Thus, we 
determine 4>j = (pjiN, 0init)- We repeat this analysis for 
different N, (/jinit, and initial configurations to obtain the 
statistics of <f>j reported in Table [J The main result of 
this analysis is that (f)j{N — > oo, (pi-nit) seems sharply de- 
fined (the fluctuations decrease in good agreement with 
the N~^^^ scaling for each value of (/jjnit, but the 
0init dependence survives a substantial increase of sys- 
tem size. Therefore, we have directly established that, for 
frictionless spheres, the jamming transition occurs along 
a continuous line rather than at a specific point in vol- 
ume fraction, as hinted in earlier simulations 1^, 22| , and 
predicted theoretically 14, lit . 

Having found several critical points for jamming, 
rather than one, we now ask whether these different 
J-points are equivalent. We have measured the static 
structure of the jammed configurations through the (par- 
tial) pair correlation functions gapir) {a, (3 = 1,2 for 
small/large particles, respectively) and the statistics of 




FIG. 2: Scaling behaviour along the continuous line of J 
points. The pressure vanishes as P ~ (0 — 0j) (left axis) 
while the number of contacts (right axis) is discontinuous, 
2 ~ Zc + V'P ~ 4' J with Zc — 6. We show data for 4 represen- 
tative points J with 4 different (pj. 



contacts between particles. Wc have also followed the 
properties of configurations approaching jamming from 
large volume fraction. As in wc find several scaling 
laws near our four J-points, for P, e and z. In Fig. [2] 
we show that for all J-points we have P ^ (cj) — (pj) and 
z = Zc + V4' ~ (f'J with Zc — Q (see Table |l|. Within our 
numerical precision, the critical properties of jamming 
transitions along the J-line are identical. 

This analysis suggests that the value of (pj is 
strongly influenced by the entire sample history (com- 
pression/decompression protocol, reorganization during 
compressions, etc.), but that the scaling properties of the 
jamming transition are much more robust. It is tempting 
to formulate the appealing hypothesis that the history 
dependence of the jammed conflgurations are contained 
solely in the volume fraction (jjj at which they jam. This 
would suggest that data analysis of the jamming transi- 
tion should drastically simplify if we eliminate (fij from 
consideration. We confirm the validity of our hypothesis 
in Fig. [T]-d, where the three sources of fluctuations dis- 
cussed above are accounted for by eliminating (f) from the 
description. By plotting e vs. z, we find that data for all 
cases nicely collapse onto a mastercurve, e ~ (z — z^)^, 
independently of the sample preparation history. Thus, 
the details of the numerical protocols are irrelevant to 
the relationship between energy and coordination num- 
ber across the jamming transition. 

Not only are scaling properties near jamming robust, 
but the special structure of the pair correlation function 
at the transition is also observed in all our samples. The 
pair correlation has a delta peak at contact, which, when 
integrated yields the average number of contacts, see Ta- 
ble HI The total number, z, and partial numbers, Zq/s, 
have the same value for all J-points, within numerical 
accuracy. Note in particular that Z12 does not evolve, 
indicating that dcmixing between species of the mixture 
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FIG. 3: (a) Statistical distribution of the number of contacts 
per particle, z, for 4 points J. (b) Power law singularity near 
contact for (/ii(r) ~ (r — cru)^ for the four points J, with 
7 = ^ shown as dashed line, (c) 311 (r) near second and third 
peaks for 4 points J. (d) As in (c) for the equilibrated fluid 
configurations at (^init. 



is not observed. In fact the full distributions of contact 
number barely change, see Fig. [3]-a. The pair correla- 
tion also has a power law singularity near contact Q , 
9ai3{r) ^ {r — CTafs)'^- This is shown in Fig. [3]-b for gii(r), 
where again the exponent is the same for the four J- 
points, 7 « 0.5. In Fig. [3]-c, we observe a subtle but 
systematic evolution with 0init of g{r) from near the first 
minimum to the third peak. Similar to the case of in- 
herent structures in systems with soft potentials, these 
subtle variations reflect the more clearly visible changes 
in the structure of the equilibrated fluid configurations 
at ^init which were used as starting points for the com- 
pressions, as shown in Fig. [3]-d. Finally, we note that 
the increase of 4>j does not result simply from a decreas- 
ing number of rattlers (particles with no contact), see 
Table H 

Our results directly establish that a unique location for 
the jamming transition cannot be expected, even when 
crystallization or demixing play no role. Rather, our find- 
ings are broadly consistent with a 'landscape' picture 
wherein the fluid explores a phase space that does not 
form the basin of a well-defined, unique ground state, 
but instead one marked by the presence of many local 
minima or metastable states, as is typical of glassy sys- 
tems. The basins predominantly sampled vary with the 
equilibrium volume fraction of the system. As in ther- 
mal systems [131 using the inherent structure formalism, 
the nature, number and evolution with volume fraction of 
these metastable states can be evaluate d |20||. a nd are de- 
scribed exactly in the mean- field limit [14l . llSj . An eval- 



uation of configurational entropy associated with these 
jammed states and a comparison with the sharp bound- 
aries predicted in mean-field theory, corresponding to the 
appearance and vanishing of a configurational entropy, 
would be valuable for this system. 

In summary, we have shown that jamming transitions 
of amorphous packings of frictionless spheres occur along 
a continuous range of volume fraction. These transition 
are sharply defined in the thermodynamic limit, and this 
finite range exists even when ordering phenomena or fric- 
tion are absent, contrasting with the idea of a unique 
jamming point. Since the location of the jamming tran- 
sition in fact results from the specific protocol used to 
study it, we see no reason why it should be reproducible 
from one experiment to another. 
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